The Width of the Confining String in Yang-Mills Theory 



F. Gliozzi", M. Pepe^ and U.-J. Wiese'= 
" Dipartimento dt Fisica Teorica, Universtta di Torino, and INFN, 
Sezione di Torino, via P. Giuria 1, 10125 Torino, Italy 
' INFN, Sezione di Milano-Bicocca, Edificio U2, Piazza della Scienza 3, 20126 Milano, Italy 
Albert Einstein Center for Fundamental Physics, Institute for Theoretical Physics, 
Bern University, Sidlerstrasse 5, 3012 Bern, Switzerland 

We investigate the transverse fluctuations of the confining string connecting two static quarks in 
(2 + l)-d SU{2) Yang-Mills theory using Monte Carlo calculations. The exponentially suppressed 
signal is extracted from the large noise by a very efficient multi-level algorithm. The resulting 
width of the string increases logarithmically with the distance between the static quark charges. 
Corrections at intermediate distances due to universal higher order terms in the effective string 
action are calculated analytically. They accurately fit the numerical data. 
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Understanding the dynamics of confining strings con- 
necting static quarks and anti-quarks in non-Abelian 
gauge theories is a great challenge in strong interaction 
physics. During its time-evolution, a confining string 
sweeps out a world-sheet whose boundary is determined 
by the world-lines of the static external quark charges. 
The string world-sheet can also be viewed as an ana- 
log of a fluctuating interface separating different phases 
of condensed matter. Just like a rough interface, the 
confining string in a non-Abelian gauge theory supports 
massless transverse fluctuations. Interestingly, the dy- 
namics of these fluctuations — known as capillary waves 
in condensed matter physics — is captured by a system- 
atic 2-dimensional low-energy effective field theory. In 
the effective theory, the string is described by a (d — 2)- 
component vector h{x,t) pointing to the location of the 
string world-sheet in the (d — 2) transverse dimensions of 
a d-dimensional space-time. Here x € [0, r] and t G [0, f3] 
are the Euclidean coordinates parameterizing the base- 
space, with r being the fixed distance between the exter- 
nal static charges and the inverse temperature /3 being 
the extent of Euclidean time. Since the string ends at 
the static quark charges, its fluctuation field obeys the 
boundary condition h{0,t) = h{r,t) = 0. The leading 
term in the action of the effective theory 

S[h] ^'^ f dt ( dx djl ■ d^h, ^ e {1, 2}, (1) 
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gives rise to a universal contribution to the static quark 
potential 

V{r)^ar-^^ + Oil/r% (2) 

where cr is the string tension. The universal Liischer term 
is proportional to the number (d— 2) of transverse direc- 
tions in which the string is fluctuating [l], Q . The effec- 
tive action from above also accounts for the string width 
Due to its transverse fluctuations, the string broad- 
ens as the quark sources are separated. More precisely 



the transverse area swept out by the flux tube increases 
logarithmically with separation. At the distance r/2 half- 
way between the external static quark sources, at leading 
order the string has the squared width 

wUr /2) log (3) 

where tq is some distance scale. This fundamental for- 
mula for the width is again universal. It applies to strings 
in confining gauge theories, to fundamental strings, as 
well as to fluctuating interfaces in condensed matter 
physics. Interestingly, until now the logarithmic behavior 
of the string width has been verifled only in the context 
of spin models. Below the critical temperature of the 
3-d Ising model, but above the corresponding roughen- 
ing transition, an interface separating the broken phases 
indeed has massless fluctuations described by the effec- 
tive theory of eq. ([T]) . Using an efficient cluster algorithm, 
the expression for the width of eq.® has been verified 
in [4|. By a duality transformation, the results obtained 
in the 3-d Ising model also apply to (2 -I- l)-d Abelian 
Z(2) lattice gauge theory, where the effective theory has 
been confirmed in more detail also using different bound- 
ary conditions [H'] . Similar results have been obtained in 
a (2 -I- l)-d Z(4) gauge theory 0. In non-Abelian gauge 
theories the measurement of the string width is computa- 
tionally very challenging; early calculations were affected 
by large statistical errors and did thus not lead to definite 
conclusions 0- 

In this paper, for the first time we verify the logarith- 
mically divergent string width for a non-Abelian gauge 
theory by simulating (2 + l)-d SU{2) lattice Yang-Mills 
theory using a very efficient multi-level algorithm Q ■ 
As we will see, the Monte Carlo data agree with the ana- 
lytic prediction at large distances r. Before we turn to the 
numerical results, we systematically work out corrections 
at intermediate distances, which arise due to higher order 
terms in the low-energy effective action. For d = 3 there 
is only one transverse dimension, and hence in this case 
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h{x,t) € R is a one-component scalar field. The lead- 
ing and next-to-leading terms in the low-energy effective 
action describing the massless transverse fluctuations of 
the string are then given by 



S[h] 



dt / dx 
Jq 
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Since the Dirichlet boundary conditions h{0,t) = 
h(r,t) = explicitly break translation invariance, one 
may have expected boundary terms to be present in the 
effective action as well j^. Remarkably, due to open- 
closed string duality, such terms are absent and the pre- 
factor of the first sub-leading term is uniquely determined 
[loj . A generalization of these arguments shows that for 
d = 3 there is only one six-derivative term which coin- 
cides with the one of the Nambu-Goto action 

The effective action of eq.(|l]) describes string fluctua- 
tions in the continuum. Before one reaches the contin- 
uum limit, the confining string in a lattice Yang-Mills 
theory is also affected by lattice artifacts. First of all, 
at very strong coupling the world-sheet swept out by the 
lattice string is rigid, i.e. it follows the discrete lattice 
steps and does not even have massless excitations. Only 
at weaker coupling, after crossing the roughening transi- 
tion, the string world-sheet supports massless excitations 
and thus becomes rough. Consequently, the effective the- 
ory is applicable only in the rough phase. 

For a string world-sheet with periodic boundary condi- 
tions in Euclidean time, the squared width of the string 
is given by 



w'^{x) = {h{x,tf 



(5) 



which is directly related to the two-point function 
{h{x,t) h{x' ,t')), with the two points {x,t) and {x',t') 
falling on top of each other. This limit leads to ultraviolet 
divergences which we regularize using the point-splitting 
procedure 

{h{x,tf) — > {h{x,t)h{x' ^ x + e,t' ^t + e')). (6) 

At the end, one sets e — >■ , e' — > and the remain- 
ing ultraviolet divergent terms are absorbed in physical 
length scales. In the leading order Gaussian approxima- 
tion of eq. , the two-point function can be conveniently 
expressed in the form 

1 . . e~"^2 -I- g"e"^2 

{h(x^t) h(x' ,t')) — — > sin n^i sin n^J^ - 
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where = nx/r, = t^x' /r, ^2 = 7r(i — t')/r with 
< ^2 < TT/S/r. The parameter q — q^t^^t depends on the 
ratio T = ij3/{2r). One then obtains the squared width 
of the string as 
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with ro = 7r|e -I- ze'l/2 and the Dedekind function 



77(r) = q 
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In the numerical simulations we put r <^ j3. Then there 
are only exponentially small corrections to the logarith- 
mic broadening of the width of the string. On the other 
hand, for r 3> /?, the inversion transformation property 



r;(r)=,7(-l/r)/^ 
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implies that 
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This shows that, at finite temperature, the confining 
string broadens linearly with the distance between the 
external static quarks 12 1. 

The non-Gaussian correction due to the first sub- 
leading term in eq.Q results from the six-point function 



(T V T 



dx / dt {dfj,hdfj_hY 



(12) 

This amplitude has two kinds of ultraviolet divergent 
terms: one proportional to log |e-|-ie'| and the other pro- 
portional to (e^ — e'^)/|e -I- «e'|^. The former is absorbed 
in the scale of the logarithm and the latter is eliminated 
by putting — e'^. Therefore the contribution to of 
the next-to-leading term does not require the introduc- 
tion of new low-energy parameters. Up to first order in 
the expansion parameter l/(f7r^), the complete expres- 
sion turns out to be 



w'{r/2) = 1 + 
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Here E2{t) is the first Eisenstein series. Just like the ac- 
tion, the operator h{x, t) also receives higher-order cor- 
rections. Remarkably, their effect vanishes in eq. (ll3|) . 
The derivation of these results will be presented in 13]. 

Lattice gauge theory is a powerful tool that allows us 
to address the question of string dynamics from first prin- 
ciples. Since it is least problematical for numerical sim- 
ulations, we consider (2 -I- l)-d SU{2) Yang-Mills theory 
on a cubic lattice of size L x L x j3, with the Euclidean 
time extent /3 determining the inverse temperature. We 
use the standard Wilson plaquette action 
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for parallel transporter variables Ux,fi G SU{2) in the 
fundamental representation of the gauge group, located 
on the links {x, fi). Here g is the bare gauge coupling and 
/t is a unit-vector pointing in the /^-direction. It should be 
noted that all physical quantities are measured in units 
of the lattice spacing which we put to 1. The partition 
function takes the form 



Z = VU e 



-S[U] _ 



n 



SU(2) 



-SlU] 



(15) 



where dUx^fi denotes the local gauge invariant Haar mea- 
sure on the link (x,ii). An external static quark located 
at the site x is represented by a Polyakov loop 



= -Tr 
2 



x+t2,2 



(16) 



a parallel transporter wrapping around the periodic Eu- 
clidean time direction. We have chosen the 2-direction 
to represent Euclidean time. When two external quark 
charges are separated along the 1-direction, the world- 
sheet of the confining string extends in the space-time 
directions /i = 1,2, while the string fluctuates in the 
transverse 3-direction. The static quark potential V{r) 
results from the Polyakov loop correlation function 



^n^r) = I- I VU <^f)^r e 



-slu] 



-PVir) 



(17) 



in the zero-temperature limit /3 — > oo. 

In order to ensure a good projection on the ground 
state of the string and in order to cover a wide range of 
spatial distances satisfying r < /3, we have simulated at 
inverse temperatures as large as /3 = 48 in lattice units. 
The spatial lattice size was L — 54 and the bare gauge 
coupling was chosen as A/g'^ — 9.0 which puts the decon- 
finement phase transition at /3c ~ 6. The numerical sim- 
ulations have been performed using the Liischer-Weisz 
technique [1, 0] with 2 levels. We measure the static po- 
tential using the Polyakov loop two-point function and 
we obtain a = 0.025897(15) for the string tension: this 
result is compatible with the value reported in 14j . 

The width of the fluctuating string is obtained from 
the connected correlation function 



C(X3) 



- (P. 



of a pair of Polyakov loops with a single plaquette 



= ^Tr[C/,,iC/, 



2; + l,2'^^+2,l^a;,2J' 



(18) 



(19) 



in the 1-2 plane, which measures the color-electric fleld 
along the 1-direction of the string as a function of the 
transverse displacement ^3. The plaquette is located at 
the site x = (r/2,0,X3) and thus measures the trans- 
verse fluctuations of the color flux tube at the maximal 




FIG. 1: The ratio {^o^rPx) / (^o^r) as a function of the 
transverse displacement X3 at fixed distance r = 19 between 
the external static quarks. The solid line is a fit of the nu- 
merical data using eg. S20\) . 



distance r/2 from the external quark charges. The corre- 
lation function at distance r = 19 is illustrated in figure[T] 
The data show the expected bell-shape of a Gaussian dis- 
tribution, but their high numerical accuracy also reveals 
small deviations. We fit the data using the ansatz 



i^O^rPx) 



l + B cxp(-a;2/T) 
1 + D exp{-xl/T) 



K, 
(20) 

where A, B, D, T, and K are fit parameters. This func- 
tion always provides an excellent flt of the data. The 
squared width of the string is then obtained as the sec- 
ond moment of the correlation function 



J dx3 xl C{x3) 
J dx3 C{X3) 



(21) 



In order to check the dependence of the measured string 
width on the ansatz of eq. ipO)) . we have also considered 
other flt functions. As long as the ansatz provided a very 
good flt of the numerical data, we have not observed a 
significant effect on the measured string width. Hence, 
systematic uncertainties due to the fitting ansatz were 
negligible compared to the statistical errors. In figure [2] 
we illustrate the dependence of the squared string width 
w'^{r/2) on the distance r between the external static 
quarks. At distances larger than r ^ 18, the Monte Carlo 
data are well represented by the effective field theory pre- 
diction of eq. (|13|) . In particular, the predicted prefactor 
l/(27rCT) of the logarithmic term log(r/ro) is confirmed by 
the numerical data. It should be noted that the string 
tension has already been determined by the flt to the 
static quark potential. The only free parameter in the flt 
of the Monte Carlo data for the width is the scale rg : we 
obtain ro = 2.26(2) such that = 0.364(3). 

It should be noted that the effective theory prediction 
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FIG. 2: r/ie squared width of the confining string uP'{r/2) 
at its midpoint as a function of the distance r between the 
external quark charges. The solid curve is a fit to the next-to- 
leading order prediction of the low-energy effective field theory 
from ea.l[T3\). 



between the external static quarks. The numerical data 
show that the normalized probability distribution is not 
significantly affected by the different probes used to ex- 
tract the color flux tube width. Early investigations in 
4-d can be found in 1151. 
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of eq. (fT3|) does not include lattice artifacts due to the 
violation of rotation invariance. Since the lattice theory 
is invariant only under discrete rotations and not un- 
der the full Poincare group, before one reaches the con- 
tinuum limit the two additional terms 2i^f^'^f^^)^ 
and J2fj.=i 2('^t^^)'^ enter the effective theory. Since these 
terms contain four derivatives, they are of next-to-leading 
order. Hence, they have no effect on the Liischer term or 
on the leading logarithmic behavior of the string width. 
The corrections due to the rotation symmetry breaking 
terms can be computed [l3| and included in the fit for the 
string width. This adds two free parameters and makes 
the fit excellent down to r ss 10. However, an approxi- 
mate estimate of the correction to the string width com- 
ing from the rotation invariant next-to-next-to-leading 
order may also significantly improve the agreement in 
the range r = 10 — 18. Since, at the moment, we cannot 
clearly disentangle these two effects, we have not taken 
into account rotation symmetry breaking corrections in 
our data analysis. 

In the numerical results discussed so far, the plaquette 
orientation was parallel to the string world-sheet. The 
plaquette acts as a probe that measures the fiuctuations 
of the confining string. The width depends on the probe 
through the value of the low-energy parameter tq ■ Choos- 
ing different probes — and thus different ways of defin- 
ing the string width — we expect, however, only small 
changes in the value of the low-energy parameter tq . In 
order to investigate this issue, we have considered dif- 
ferent orientations of the plaquette in the definition of 
eq. ()18|) . In figure [3] we show the normalized probabil- 
ity distribution C{xy,)/ J dx^, C{x3) for the 3 different 
orientations of the plaquette at fixed distance r 12 



FIG. 3: Probability distribution C{x3)/ J dxg C(xs) using the 
3 possible orientations of the plaquette in eg. \18\) . The data 
points for the 2 orientations orthogonal to the string world- 
sheet are put at the midpoint between 2 lattice sites. 
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